Derealization of ultracold atoms in a disordered potential due to light scattering 
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We numerically study the expansion dynamics of ultracold atoms in a one-dimensional disordered 
potential in the presence of a weak position measurement of the atoms. We specifically consider 
this position measurement to be realized by a combination of an external laser and a periodic 
array of optical microcavities along a waveguide. The position information is acquired through the 
scattering of a near-resonant laser photon into a specific eigenmode of one of the cavities. The 
time evolution of the atomic density in the presence of this light scattering mechanism is described 
within a Lindblad master equation approach, which is numerically implemented using the Monte 
Carlo wave function technique. We find that an arbitrarily weak rate of photon emission leads to a 
breakdown of Anderson localization of the atoms. 



I. INTRODUCTION 

The realization of potentials with controlled disorder 
for ultracold atoms by means of optical speckle fields 
P, or Dichromatic optical lattices [|| has recently led 
to the observation of Anderson localization with Bosc- 
Einstein condensates [3, H[. In those experiments, atomic 
Bose-Einstein condensates, prepared in a harmonic trap, 
were released into one-dimensional optical waveguides 
which were superimposed with disordered potentials real- 
ized with speckle fields [|| as well as with bichromatic op- 
tical lattices Absorption images of the atomic cloud 
after the expansion process within the waveguide clearly 
revealed an exponential decrease of the average atomic 
density with the distance from the center of the former 
trap, which is the characteristic signature of Anderson 
localization Q. While interaction effects did not play a 
role in those pioneering experiments, more recent stud- 
ies specifically focus on the interplay of atom-atom in- 
teraction and localization in disordered potentials (e.g. 
Ref. 0,11). Current research directions include the ex- 
ploration of Anderson localization with ultracold atoms 
in three spatial dimensions with the particular aim 
to study the Anderson metal- insulator transition (Tfjj . 

Clearly, a key condition for the observability of Ander- 
son localization with ultracold atomic gases is the overall 
coherence of the atomic cloud. Any mechanism of deco- 
herence would compromise the phenomenon of destruc- 
tive wave interference that lies at the heart of Anderson 
localization Q and thereby give rise to derealization. 
This also concerns any in-situ monitoring of the evolu- 
tion of the atomic cloud during its expansion, by interme- 
diate measurements of the positions of atoms. Evidently, 
the strong refocusing of the atomic wavefunction that re- 
sults from a precise position measurement would destroy 
the coherence of the atom, enhance its kinetic energy, 
and eventually let the atom behave as a classical particle 
when being performed several times. 



■ \ 


1=1 




1=1 










1 " \ 




















^= 









FIG. 1: (Color online) Sketch of the configuration under 
consideration. An atom, which is initially localized within 
a harmonic trap at the center, is expanding within a one- 
dimensional waveguide (indicated by the green horizontal 
lines) to which a disordered potential is superimposed. A 
periodic array of optical microcavities (symbolized by the red 
arcs) is used to measure the position of the atom, on a length 
scale that is comparable to its localization length (indicated 
by the shadow plot in the waveguide, which shows the den- 
sity of a localized state). For this purpose, the waveguide is 
considered to be irradiated by a near-resonant laser beam, 
which may induce spontaneous emissions of photons into one 
of the cavities, possibly to be measured by photodetectors. 
The sketch is to scale with the parameters considered in this 
study, as far as the horizontal length scales are concerned. 



The situation is less obvious for "weak" measurement 
processes, in which the position of the atom is deter- 
mined with a large spatial uncertainty that is of the order 
of the expected localization length within the disordered 
potential. Such weak measurements might still preserve 
coherence to a certain extent, while, at the same time, 
providing some rough in-situ information on the position 
of the atom. One could, for this purpose, conceive e.g. 
a periodic array of optical microcavities placed around 
the waveguide in which the atoms propagate, as depicted 
in Fig. Q] A near-resonant laser beam which irradiates 
this configuration can be used to transfer the atoms to 
an electronically excited state, from which they can re- 
lax to the ground state by a spontaneous emission of a 
photon into one of the cavities, which in turn could be 
measured by photodetectors placed behind the cavities. 
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The whole configuration could possibly be fabricated on 
"atom chips" in which case the disordered poten- 

tial could arise from imperfections in the current-carrying 
wires that generate the magnetic waveguide potential of 
the atoms fl2l Il3j. Our results below are, however, more 
general and we expect to see the same effects on the lo- 
calization in the presence of any mechanism of similar 
position measurement. 

The aim of this study is to investigate to which ex- 
tent this approximate realization of a "Heisenberg micro- 
scope" gives rise to derealization of an atom in a one- 
dimensional disordered potential. For the sake of sim- 
plicity, we shall restrict our consideration to the propa- 
gation of one single atom, and thereby discard collective 
processes arising within Bose-Einstein condensates due 
to atom-atom interactions or superradiance. We shall, 
moreover, assume that the atom will emit photons into 
one single mode of the cavities only. Such an emission 
will then give rise to a recoil that is mainly perpendicu- 
lar to the direction of propagation of the atom and does 
therefore not dramatically enhance its longitudinal ki- 
netic energy. We neglect effects of transverse excitations 
within the waveguide due to this recoil and assume that 
neither the effective waveguide confinement nor the dis- 
ordered potential are affected by temporary populations 
of the excited electronic state of the atom. 

The dynamics of the atom is modeled via a one- 
dimensional master equation for its density matrix p(t), 
which can be unraveled using the Monte-Carlo wavefunc- 
tion technique flil Il5j . This master equation accounts 
both for the coherent motion within the disordered po- 
tential and the incoherent scattering of photons [l6| . We 
shall, in Section [TH first account on the expansion and lo- 
calization dynamics of a single atom in a one-dimensional 
disordered potential in the absence of any decoherence 
mechanism. In Section IIII| we outline the Monte-Carlo 
wavefunction approach that is used to integrate the mas- 
ter equation for the special case of an atom that propa- 
gates in a homogeneous, disorder-free waveguide. Deco- 
herence and disorder are finally put together in Section 
IIV1 in which we discuss the expansion of an atomic wave 
packet in the presence of disorder and spontaneous emis- 
sion. We show that even very rare position measurements 
of the atom give, on average, rise to a gradual dereal- 
ization of the wave packet, and we provide numerical 
evidence for superballistic expansion in the presence of 
strong emission rates. 



II. WAVE PACKET EXPANSION IN DISORDER 

In this section, the expansion of an initially trapped 
wave packet in a weak one-dimensional disordered poten- 
tial is discussed. For the sake of simplicity, we model the 
disorder by a Gaussian correlated random potential V{x) 
defined along the a:-axis, with the properties V(x) = 
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FIG. 2: (Color online) Time evolution of the expansion pro- 
cess of a Gaussian wavefunction, averaged over 100 disorder 
realizations, for the strength U — 0.1 hu and the correlation 
length a = 0.2 ao of the disordered potential (the mean ini- 
tial kinetic energy of the wavefunction is E = 0.25 hu, where 
ao and u are the oscillator length and the frequency of the 
harmonic confinement potential, respectively) . Shown are the 
disorder-averaged probability densities p(x) for cut = 20, 200, 
400, 600, and 800 (with the trap opening at t = 0). The solid 
line displays the analytical prediction ((SJ) which is found to 
be in good agreement with the numerical density distribution 
at large distances \x\ S> ao, apart from a global prefactor of 
the order of 2. 

and 

V(x)V(x') = U 2 exp[-(z - x') 2 /(2a 2 )} (1) 

for the mean spatial correlation function. Here, U charac- 
terizes the typical size of the fluctuations of the potential, 
and the correlation length a controls the average width 
of fluctuations. 

In Fig. [2] we show the time evolution of the disorder- 
averaged spatial density of wave packets propagating in 
such disorder configurations. These wave packets are ini- 
tially prepared in the ground state of a harmonic trap 
with the oscillator length ao = \J h/muj. After the trap- 
ping potential is switched off, the wave packet expands 
within the disordered potential until it approaches, on 
average, a stationary profile. The convergence to the 
average density distribution happens faster at the center 
than in the wings. This is a consequence of the quadratic 
growth of the localization length as a function of the wave 
vector, as described in Eq. |0J) below. 

The final density profile shown in Fig. [2] is fairly well 
reproduced by a theory as described, for example, in 
Ref. (l7j . which is based on the assumption that the 
asymptotic probability distribution is an incoherent sum 
of individually localized plane waves with momentum p. 
This consideration yields the spatial density 

Pioo(x) = J dpgljexpH*!/^)] (2) 

where po denotes the momentum density of the wave 
packet at the initial time t = 0. The key ingredient for 
the evaluation of Eq. ^ is the localization length 
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that can be calculated using diagrammatic theory [18| as 
£(p) = 21 nip), with Ib the Boltzmann mean free path 
(see Refs. [US, [l9| for other approaches). For the Gaus- 
sian correlated random potential under consideration, we 
obtain 
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In the regime of short correlation lengths a <C h/p, we 
can approximate exp[2(per//i) 2 ] ~ 1 and the localization 
length depends only on the effective strength U 2 a of the 
disorder. Using 



Po(p) 



exp[-(a p/h) 2 



(4) 



and introducing the characteristic localization length 
scale of the wave packet as £o = £,{h/a.o), we then ob- 
tain the prediction 



P\oc{x) 
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: CXp 
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for the localized density. As shown in Fig. this ap- 
proximate expression is, apart from a global prefactor, 
in good agreement with the numerically computed mean 
density at the final time t = 800/w. 

In the above numerical simulations, we effectively as- 
sumed that the atomic cloud is prepared in a clean har- 
monic trap in absence of any disorder. At t = the 
trapping potential is suddenly switched off and the dis- 
order is ramped on at the same time. The initial state is 
then a perfect Gaussian wavefunction [see Eq. ([4])] which 
expands within the disordered potential. This procedure 
is, in general, not precisely in accordance with expansion 
experiments on Anderson localization such as Ref. Q in 
which the disordered potential is already present during 
the formation of the Bosc-Einstcin condensate in the har- 
monic trap. The initial state of the atomic wavefunction 
is, in that case, given by the ground state of an effec- 
tive trapping potential that consists of a harmonic con- 
finement modulated by the disorder. A numerical com- 
parison of these two expansion scenarios, however, dis- 
plays no significant difference in the asymptotic density 
profile for the case of weak disordered potentials with 
U ~ 0.1 hui and a = 0.2 ag. 

A convenient numerical observable for measuring lo- 
calization is the participation ratio [2(| which for a wave 
packet with the density p(x, t) is defined by 



Pr(i) 



dx[p(x, t)] 



(6) 



In practice, Pr(£) represents a measure for the spa- 
tial extent of the wave packet, yielding large values for 
rather extended distributions p(x, t) and going to zero 
for strongly peaked wavcf unctions. It therefore exhibits 
a similar behavior to the spatial root mean square (rms) 
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FIG. 3: (Color online) Root mean square (rms) width Aa; = 
■yj (x 2 ) — (x) 2 (red line) and participation ratio Pr(t) (black 
line) as a function of the evolution time t, showing the ex- 
pansion and localization of a wave packet for the disorder 
strength U = 0.15 huj and the correlation length a — 0.2 ao. 
The dashed lines show, for comparison, the rms width and 
the participation ratio of a free wave packet that expands in 
the absence of disorder. 



width Ax = yj (x 2 ) — (x) 2 of the wave packet. This lat- 
ter quantity is, however, rather sensitive to the evolution 
of the (experimentally inaccessible) wings of the wave 
packet. This is shown in Fig. [3] where we display the 
time dependence of the disorder-averaged rms width and 
participation ratio. While the rms width continuously in- 
creases with time, due to the long-time dynamics in the 
wings of the averaged density distribution (see Fig. [2]), 
the participation ratio, which is much less sensitive to 
the behavior of the wings, saturates at a finite length 
scale. This length scale can be used in order to define an 
effective localization length L wp of the wave packet. 

Fig. [4] shows the time evolution of the participation 
ratio for different initial kinetic energies E, resulting from 
different confinement frequencies of the initial trapping 
potential. A linear increase of the participation ratio 
with E, corresponding to a linear increase of the effective 
localization length L wp of the wave packet, is found for 

Eq < E < 2Eq. 



III. MASTER EQUATION DYNAMICS 

To account for spontaneous emissions of photons into 
the cavities, we model the dynamics of the atom via a 
one-dimensional master equation for its density matrix 
p(t), including coherent interactions with a disordered 
potential and the incoherent scattering of light [16[ . This 
master equation is given by 



d , , 



- h [Hr P {t)] 



7c ff r ^(c q p(t)ct-p(t)) 



(7) 

Here, H = # h V(x) describes the Hamiltonian for 

a particle that propagates in the disordered potential. 

is the decay or jump operator representing 
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one spontaneous emission event, which exerts a recoil 
on the atom with longitudinal momentum hq which is 
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FIG. 4: (Color online) Left panel: participation ratio as a 
function of time for different initial kinetic energies E, with 
Eo — 0.25 fiw. The disordered potential is characterized by 
the parameters U = 0.15 hjJ and a = 0.2 do- The right panel 
shows the average of the participation ratio within 400 < uit < 
1000 as a function of the energy E. This average defines the 
asymptotic localization length L wp of the wave packet. L wp 
is found to increase approximately linearly with the initial 
kinetic energy. 



assumed to be equidistributed between —hk and +hk. 
This model considers off-resonant inelastic scattering, in 
which a laser couples the electronic ground state to an ex- 
cited state from which spontaneous emission back to the 
ground state can occur. It assumes a low spontaneous 
decay rate 7 as compared to the detuning 5 of the laser 
with respect to the intra-atomic transition frequency, and 
a low Rabi frequency v for laser-induced transitions be- 
tween the ground state and the excited state as com- 
pared to the spontaneous decay rate 7, i.e. we assume 
v -C 7 <C 6. We then obtain 7 e g = 7^ 2 /(7 2 + 4(5 2 ) as the 
effective decay rate that enters the master equation ([7| 

MM. 



To solve the time evolution generated by the mas- 
ter equation, we employ the Monte Carlo wave func- 
tion method [3, [Hi- Here, the evolution of the den- 
sity matrix is decomposed into the non-unitary evolu- 
tion of a large number N = 100 of wave functions. A 
single trajectory i = 1...N evolves according to 

ihd t \ip) t = H c ff\ijj)i with H eS = H - ijes/2, until the 
exponentially decaying norm ||V')i| 2 = e~ 7off * equals a 
random number chosen between and 1. At this point, 
a jump operator C q acts on the Monte Carlo wave func- 
tion: \ip(t + St))i = C q \ip(t))i. This jump operator is 
determined by randomly choosing q from the interval be- 
tween —k and k. 

To relate this light scattering process to the position 
measurement under consideration, we note that the Lind- 
blad master equation is invariant under unitary trans- 
formations on the set of decay operators. Indeed, it 
was shown in Ref. [23j that the Fourier transformation 
J_ 1 du exp{iukv\/2) C u with integers v € Z allows one 




FIG. 5: Center of mass position (x) (solid line) for a single 
trajectory of an atom that subject to light scattering with 
k — 0.07/ao. The dashed lines indicate the root mean square 
(rms) width Ax = -*/ (x 2 ) — (x) 2 , i.e. they correspond to the 
lines (a;) ± Ax. Photon emissions occur at u)t = 60, 79, 113, 
and 148. 



to switch to decay operators 

-sin(fci — -|) 



C V = V2- 



kx 



(8) 



In this picture, the application of the decay operator in- 
duces a localization of the wavefunction within a spa- 
tial region whose extent is of the order of k" 1 . For the 
sake of simplicity, we assume that these decay operators 
exactly correspond to the longitudinal structure of the 
cavity modes into which the atom may emit the photon. 
The spatial period of the array of cavities is then given 
by A = 27r/fc. 

In Fig. [5l we show the time evolutions of the ex- 
pectation value of the position (x) and its rms width 
Ax = y (^ 2 ) ~ ( x )' 2 for a single quantum trajectory in 
free space, V(x) = 0. In this particular trajectory, the 
first spontaneous emission took place in one of the two 
wings of the wave packet, which is mainly constituted by 
plane-wave components with high momenta. The subse- 
quent localization process projects the wavefunction on 
those high-momentum components, which gives rise to a 
permanent drift. The rms width, however, remains small 
during this evolution, which is due to the fact that the 
atom emits photons at a rate that is faster than the in- 
verse dispersion time of the wave packet. The rms width 
would freeze for sufficiently high emission rates, which is 
reminiscent of the quantum Zcno effect. 

It is of great advantage to work in a regime where k 
is small compared to fluctuations of the density matrix 
in momentum space p(p,p',t). To study the momentum 
density distribution, we can then approximate the inte- 
grand of Eq. by its Taylor expansion to first order, as 
done in Ref. [24[ . Taking the integral over k leads to the 
diffusion equation 
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(9) 



for the diagonal elements of the density matrix, with the 
effective diffusion coefficient D = j c gk 2 /6. Hence, the 
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FIG. 6: (Color online) Master equation dynamics of the par- 
ticipation ratio for different effective emission rates 7 e ff. The 
longitudinal momentum of the emitted photon is k — 0.07/ao 
in the left panel and k — 0.035/ao in the right panel. The 
linear growth of the participation ratio for cut > 100 reflects 
the growth of the mean kinetic energy due to photon scat- 
tering events, which is found to increase from E — 0.25 huj 
at t = to E ~ 0.5 fej at ujt = 1000 for k = 0.07/a and 
7cff = 0.05oj (blue curve in the left panel). The disordered 
potential is characterized by the strength U = 0.15 ftio and 
the correlation length a = 0.2 ao. 



wave packet will undergo diffusive spreading in momen- 
tum space. Noting that the variance of the momentum 
distribution is nothing but the kinetic energy, we obtain 



(t) 
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h 2 k 
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(10) 



for the growth of the mean kinetic energy of the wave 
packet. 



IV. DISSIPATIVE EXPANSION IN DISORDER 

Having introduced the necessary tools, we now study 
wave packet expansion in the presence of disorder and 
dissipation. In Fig. [6] we plot the participation ratio as 
a function of time for different effective emission rates 
Teg. In accordance with the sketch shown in Fig. [TJ we 
have chosen the photon wavelength to be very long com- 
pared to the initial extension ao of the wave packet. The 
amount of kinetic energy given to the wave packet at each 
emission event is thereby rather reduced. 

The most important observation is a derealization of 
the wave packet at any emission rate. Instead of satu- 
rating to a stationary value, the participation ratio lin- 
early increases with time after the typical time scale that 
is needed for developing an Anderson-localized density 
profile in the absence of spontaneous emission. Quanti- 
tatively, this linear growth is very different from a bal- 
listic expansion in free space, which takes place with 
much faster expansion velocities (see the dashed lines in 
Fig. It is also different from simple diffusion which 



one would naively expect to prevail for a quantum parti- 
cle that propagates within a disordered potential in the 
presence of a decoherence mechanism. We attribute this 
difference to the fact that the spontaneous emission of a 
photon gives rise to a recoil of the atom and thereby in- 
creases its energy. Hence, the effective diffusion constant 
should also gradually increase with time. 

It is, in this context, interesting to note that the ex- 
pansion velocity dPi/dt depends only on the product of 
the effective rate of emission j c g and the recoil energy 
h 2 k 2 1 (2m). This can be seen by comparing the two blue 
lines in the two panels of Fig. [S] showing expanding par- 
ticipation rates for k = 0.07/a and 7 e ff = 0.05w (left 
panel) as well as for k = 0.035/ao and 7eff = 0.2a;. There 
appears, furthermore, no change in the behavior when we 
tune the rate of emissions across the scale l/7i oc , with 
T\ oc the time at which the unperturbed evolution shows 
localization. 

It is tempting to relate the linear increase of the par- 
ticipation rate with time to the combination of a linear 
growth of the kinetic energy due to spontaneous emission 
with the approximately linear scaling of the wave packet's 
localization length with its mean kinetic energy in the ab- 
sence of spontaneous emission, as shown in Fig. [4] This 
reasoning essentially assumes that in between two subse- 
quent spontaneous emission events the wave packet has 
enough time to approach its asymptotic stationary profile 
within the disordered potential. Extracting from Fig. [4] 
the approximate scaling Pr/an ~ IOOE/Eq and using 
dE/dt = h 2 k 2 j e g/(6m) for the growth rate of the energy 
according to Eq. ([TTJ]) . we obtain the prediction 



dPr 
~dt 
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ao dE 
E~n ~dt 
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-y c fiH 2 k 2 a 
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(11) 



for the expansion velocity dPv/dt of the participation 
rate, using E = 0.25huj. 

Figure [7] shows, however, that this expansion velocity 
increases more strongly with the rate of increase of the 
kinetic energy than predicted by Eq. (fTTj) . As a matter 
of fact, dPr/dt is found to scale as a square root of dE/dt 
in the parameter regime in which we carried out our nu- 
merical investigations. One may attribute this behavior 
to the fact that the above reasoning rather applies to 
an individual quantum trajectory in the spirit of Fig. [5] 
The energy of the wavepacket corresponding to each indi- 
vidual trajectory increases linearly and its participation 
ratio increases on average as described by Eq. (fTTj) . How- 
ever, while different trajectories describe similar narrow 
wavepackets, each wavepacket will be centered around a 
different point in space. Thus the full (incoherent) den- 
sity will be spreading faster over a larger region than a 
single wavepacket (as is obvious from Fig. [5] for the case 
of disorder- free propagation). This effect is obviously not 
accounted for in the considerations leading to Eq. (fTTj) . 

Let us finally investigate the regime of strong dis- 
sipation for which it is expected that the expansion 
becomes independent of the disordered potential. In 
Fig. [5] the momentum recoil is set to k = 0.35/ao, 
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FIG. 7: (Color online) Expansion velocity dPr/dt of the 
participation ratio as a function of the energy growth rate 
dE/dt = ~f eS h 2 k 2 /{6m). The data are extracted from Fig. \E\ 
as well as from other calculations using different parameters 
for 7 c ff and k, through linear regression of the participation 
ratio within 100 < cot < 1000. As confirmed in the log-log 
plot shown in the inset, dPr/dt scales as the square root of 
dE/dt: we have dPr/dt ~ a(dE /dt) 1/2 a /h 1/2 with the fit- 
ted proportionality constant a ~ 15.5, as indicated by the 
dashed line. The dash-dotted straight line in the main panel 
represents the prediction of Eq. (|11|) . 
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FIG. 8: (Color online) Master equation dynamics of the rms 
width Ax at strong recoil k — 0.35/ao for different effective 
emission rates y e H- The disordered potential is characterized 
by the parameters U — 0.15 hw and a — 0.2 ao- We note that 
the growth of Ax for y e B = 0.05 w (green curve) coincides 
with the one of a free ballistic expansion, while superballistic 
expansion is encountered for larger emission rates. 



drive the system beyond the free ballistic growth into a 
superballistic regime, which is the valid limit of a driven 
expansion in free space [24fj . As can be seen in Fig. [HI 
this superballistic regime sets in beyond j e s = 0.05 a; for 
k = 0.35/ao. 

V. CONCLUSION 

In summary, we have shown that even a very weak 
rate of photon scattering gives rise to a breakdown of 
Anderson localization of an atom that propagates in a 
one-dimensional disordered potential. This breakdown 
is most conveniently quantified in terms of the disorder- 
averaged participation ratio of the atomic density, which 
represents a measure for the spatial width of the atomic 
wave packet. While this participation ratio saturates, 
within a characteristic time scale, to a finite value in the 
case of a perfectly coherent expansion process within the 
disordered potential, it is found to linearly grow with 
time beyond that time scale in the presence of sponta- 
neous photon scattering. 

This growth behavior imposes strong limits for the ob- 
servability of Anderson localization in the presence of a 
weak position measurement of the atom. However, an ex- 
perimental realization of a "Heisenberg microscope" for 
cold atoms according to the scheme displayed in Fig. Q] 
might nevertheless be of interest as it allows one to study 
in more detail the interplay of disorder and measurement- 
induced delocalization phenomena not only for a single 
atom, but also (and this more naturally) for a Bosc- 
Einstcin condensate in which the atoms interact with 
each other. For this purpose, an integrated setup on 
atom chips appears as the most convenient realization of 
such a Heisenberg microscope for atomic gases. 

Finally, we expect similar findings in the presence of 
other mechanisms that can behave as a position measure- 
ment of the propagating atom. Such mechanisms include 
noise on the lattice beams as well as collisions with back- 
ground gas atoms, to mention two examples. Undesired 
effects of this type are therefore also expected to induce 
a delocalization of the atom in the disorder potential. 
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